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Quadrupolar effect and rattling motion 
in heavy fermion superconductor PrOs4Sbi2 
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The elastic properties of a filled skutterudite PrOs4Sbi2 with a heavy Fermion superconductivity 
at Tc = 1.85 K have been investigated. The elastic softening of (Cn — Ci2)/2 and C44 with lowering 
temperature down to Tq indicates that the quadrupolar fluctuation due to the CEF state plays a 
role for the Cooper paring in superconducting phase of PrOs4Sbi2. A Debye-type dispersion in 
the elastic constants around 30 K revealed a thermally activated rattling due to the off-center 
Pr-atom motion obeying r — Toexp{E/kBT) with an attempt time tq = 8.8 x 10~^^ sec and an 
activation energy E = 168 K. It is remarkable that the charge fluctuation of the off-center motion 
with r23 symmetry may mix with the quadrupolar fluctuation and enhance the elastic softening of 
(Cii - Ci2)/2 just above Tc. 

Keywords: PrOs4Sbi2, heavy fermion superconductor, quadrupolar effect, rattling motion 
PACS numbers: 74.70.Tx, 71.27.-|-a, 62.20.Dc 



o 



> 
in 

o 
o 



o 



X 



The rare-earth cubic compounds based on Pr'^+ ions 
have received much attention because various unusual 
properties are expected at low temperatures. The sys- 
tem with a non-Kramers Fa doublet possessing two 
quadrupoles 0§ = (2 ~J^~ Jy )/V3 and Oj^ 
favors a quadrupole ordering. We refer the Fa ground 
state systems as a metallic compound PrPbs showing 
the antiferro-quadrupole ordering at Tq = 0.4 and a 
semiconductor PrPtBi undergoing fcrro-quadrupole or- 
dering at Tq = 1.2 Km. PrSb is known as a singlet 
ground state system^. The elastic constant (Cn — Ci2)/2 
is responsible for the quadrupolar susceptibility of O2 
and O2 with T3 symmetry, while C44 is for the sus- 
ceptibility of O 



JyJz JzJyj Ozx — 'JzJx JxJz^ 

Oxy = JxJy + JyJx with Fs Symmetry. The softening 
of (Cii — Ci2)/2 and C44 is a useful prove to clarify the 
quadrupolar effects of Pr-based compounds. 

Recently, Bauer et al. have found a new-type of the 
heavy Fermion superconductor in a filled skutterudite 
PrOs4Sbi2 with space group T^ (Im3)^. The heavy 
Fermion state with a large specific heat coefficient 7 = 
750 mJ/mol-K^ of PrOs4Sbi2 exhibits the superconduct- 
ing transition at Tc = 1.85 K associated with a large 
jump AC/Tp ~ 500 mJ/mol-K^. A sign of the dou- 
ble transition in the specific heat has been found^. The 
thermal transport measurement in fields suggests the two 
distinct superconducting phases in PrOs4Sbi2^. The nu- 
clear spin relaxation rate 1/Ti of Sb indicates uncon- 
ventional superconductivity possessing neither a coher- 
ence peak nor a T'^-power law^,. The muon spin re- 
laxation in PrOs4Sbi2 yields a penetration depth in- 
dicating a new-type of energy gap^. The odd-parity 
Cooper pairing mediated by the quadrupole fluctuation 
is argued as unconventional heavy Fermion supercon- 



ductivity in PrOs4Sbi2^. Because the magnetic sus- 
ceptibility is rather silent to distinguish non-magnetic 
F23 doublet from Fi singlet, it has not been settled 
whether the CEF ground state of PrOs4Sbi2 is F23 dou- 
blet or Fi singlet^i^. The measurement of (Cn — Ci2)/2 
and C44 responsible for the quadrpolar susceptibility in 
PrOs4Sbi2 is a central issue to clarify the CEF state and 
the interplay of the quadrupolar fluctuation to the super- 
conductivity in PrOs4Sbi2. 

The reduction of the thermal conductivity in flUed 
skutterudites RM4Sbi2 (R: La or Ce. M: Fe or Co) is 
caused by a rattling motion due to a weakly bounded 
rare-earth ion in an oversized cage of Sb-icosahedronii. 
The filled skutterudites with the cage are favorable for 
the thermoelectric device possessing a high coefficient of 
merilii^. The ultrasonic measurements are generally use- 
ful to observe the rattling motion or off-center tunneling 
motion. We refer the rattling motion in clathrate materi- 
als SrgGaieGeaG^^ and Ce3Pd2oGe6^'*, and an off-center 
tunneling of OH-ion doped in NaCli^i^^. Recently, our 
experiment on C44 in La3Pd2oGe6 revealed the ultrasonic 
dispersion around 20 K due to the rattling motion and 
elastic softening below 3 K due to off-center tunneling 
motion in cageii. It became to be sure that the rat- 
tling and tunneling are common features of the clathrate 
compound with oversized cage. Quite recently a small 
off-center displacement of Pr-ion in PrOs4Sbi2 has been 
observed by X-ray absorption measurementsi^. The rat- 
tling motion in the clathrate compound PrOs4Sbi2 with 
the heavy Fermion superconductivity has not been clar- 
ifled yet. 

The single crystal of PrOs4Sbi2 with a length of 1.2 
mm along the [110] direction for the present ultrasonic 
measurements was grown by a flux method. The ul- 
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FIG. 1: Temperature dependence of (Cn — Ci2)/2 in 
PrOs4Sbi2 measured by ultrasonic wave of 17 MHz. The 
anomaly around 30 K originates from the rattling due to 
the Pr-ion off-center motion. The softening of (Cn — Ci2)/2 
below 20 K down to superconducting point Tc = 1.85 K is 
due to the quadrupolar fluctuation of the CEF states. The 
solid line and dashed line are fits by the quadrupolar suscepti- 

(2) (2) 

bility xq for r23-r4 and ri-r4 models, respectively. Inset 
shows the detail around Tc- 



trasonic velocity v was detected by a phase compara- 
tor based on a mixer technology. The piezoelectric 
LiNbOs transducers of i-cut and 36°^^cut were used for 
the measurements of transverse and longitudinal ultra- 
sonic waves, respectively. The elastic constant C — pv'^ 
of PrOs4Sbi2 with a lattice constant a = 0.930311 nm 
was estimated by the density p — 9.75 g/cm"^. A '^He- 
evaporation fridge down to 500 mK was employed. 

Figure 1 shows temperature dependence of (Cn — 
Ci2)/2 of PrOs4Sbi2, which was obtained by the trans- 
verse wave of 17 MHz propagating along k = [110] with 
polarization u = [110]. This (Cn — Ci2)/2 mode is as- 
sociated with the elastic strain e„ = Sxx — £yy The 
increase of (Cn — C\i)l2 around 30 K of Fig. 1 origi- 
nates from the Debye-type dispersion, where the ultra- 
sonic wave frequency w coincides with a relaxation time r 
of the system as lot = 1. A relatively large lattice param- 
eter a=0. 930311 nm in PrOs4Sbi2 may lead to the rat- 
tling motion of an off-center Pr-ion in an oversized cage 
of Sb-icosahedron. The resonant scattering of the ultra- 
sonic wave by the rattling motion of Pr-ion over a po- 
tential hill brings about the Debye-type dispersion. The 
experimental determination of a relaxation time r of the 
rattling is discussed latter. 

A remarkable softening of (Cn — Ci2)/2 below 20 K 
in Fig. 1 has been found with decreasing temperature. 
As shown in inset of Fig. 1, the softening of (Cn — 
C'i2)/2 turns up around the superconducting transition 
Tc = 1.85 K. The quadrupole-strain interaction, -ffqs = 



~ Si .9rOr7(«)er7, and the inter-site quadrupole in- 
teraction, i?QQ — — 9r{OT-y)Or-y{i), give rise to the 
elastic softening as Cp = C" — A^3rXQ/(l — ffrXo)^- 
Here means a sum over rare earth sites in unit vol- 
ume. XQ is a quadrupolar susceptibility consisting of the 
Curie-term for diagonal parts and the Van Vleck-term 
for off-diagonal ones. The coupling of the quadrupole O2 
with r23 symmetry to the elastic strain Sy is relevant for 
the softening in (Cn - Ci2)/2 below 20 K of PrOs4Sbi2. 

The CEF Hamiltonian for the Pr'^+ ion with the site 
symmetry Th is written as Hcef — 3^0^+3^0^+310^, 
where Oi ^ O^ + 50|, Oe = - 2101 and = 

Ol - Ol'^^. Two different types of the CEF models of 

(2) (2'] 
r23-r4 and ri-r4 for PrOs4Sbi2 have been proposed 

so fari^iSi. The solid line in Fig. 1 based on the doublet 
model of r23(0 K), rf^(8.2 K), ri^^(133 K), and ri(320 
K) for 3l = 6.75 x lO'^ K, 3^ = -1.23 x 10-^ K, and 
3l = -0.12 X 10^2 K ^ith jgra:,! = 97 K and an inter- 
site coupling gp^^ = —0.27 K reproduces the softening of 
(Cii — Ci2)/2 mostly proportional to reciprocal temper- 
ature. The doublet model seems to be favorable for the 
softening of (Cn — Ci2)/2. 

Recently, Kohgi et al. have proposed a singlet ground 

state CEF model of ri(0 K), ri^^(7.9 K), ri^^(135 K), 

and r23(205 K) with Bl = 2.37 x lO'^ K, = 1.32 x 

10-3 K, 3l = 1.08 X 10-2 The dashed line for 

IsrasI = 79 K and gp^^ = 0.22 K based on the singlet 

model also reproduces the softening of (Cn — Ci2)/2 

except for a small deviation below a minimum around 

3.5 K in the fitting. The off-center motion Pr ion 

with r23-symmetry may enhance the elastic softening of 

(Cn — Ci2)/2 just above Tc, that will considerably renor- 

malize the one-ion susceptibility of the dashed line in Fig. 

(2) 

1. Even in the case of the ri-r4 model, the charge fluc- 
tuation due to the off-center motion may reproduce the 
softening of (Cn — Ci2)/2 just above Tc proportional 

to reciprocal temperature. In order to settle the alter- 

(2) (2) 

native CEF model of r23-r^ ' or Fi-F^ ^ further experi- 
ments are necessary. The Debye-type dispersion was em- 
ployed to reproduce the anomaly around 20 K of solid 
and dashed lines in Fig. 1. We discuss this point again 
in Fig. 3. 

In Fig. 2, we show temperature dependence of C44 
obtained by the transverse wave propagating along [110] 
with polarization along [001]. The softening of C44 be- 
low 60 K is described in terms of the quadrupolar sus- 
(2) 

ceptibility for the F^ -type quadrupole. The solid line in 

(2) 

Fig. 2 is responsible for the r23-r4 model with param- 
eters |5r4l = 34.4 K and g'j.^ = -0.002 K. The dashed 

(2) 

line is a fit for the ri-r4 ' model with \gr^ | = 70 K and 
,gp^ — —0.07 K. Because the quadrupolar susceptibihty of 
C44 for both models is dominated by the Van Vleck term 
responsible for off-diagonal processes, the determination 
of the CEF state by C44 is rather indirect as similar as 
the magnetic susceptibility. It should be noted that no 
sign of the ultrasonic dispersion has been found in C44 
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FIG. 2: Temperature dependence of the elastic constant C44 
of PrOs4Sbi2 measured by ultrasonic wave of 74 MHz. The 
softening of below 60 K down to superconducting point Tc = 
1.85 K is described in terms of the quadrupolar susceptibility 
Xq of solid line for r23-r4 model and dashed line for ri-r4 
model. Inset shows the detail around Tq- 
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FIG. 3: Temperature dependence of Cl = (Cn + C12 + 
2C44)/2 in PrOs4Sbi2 measured by ultrasonic waves of 17.8, 
52.0 and 87.7 MHz. Arrows indicate the temperatures that 
the relaxation time r of Pr-ion rattling coincides with the 
sound wave frequencies w as ixJt = 1. Inset shows tempera- 
ture dependence of relaxation time. 



around 30 K. 

Figure 3 represents Cl = (Cn + C12 + 2C44)/2 of 
PrOs4Sbi2 obtained by the longitudinal wave along the 
[110] direction. A remarkable frequency dependence of 
17.8 MHz, 52.0 MHz and 87.7 MHz in Fig. 3 is described 
in terms of the Debye-type dispersion as, C\^{lo) = 
Cl(oo) - {Cl(oo) - Cl(0)}/(1 + Here w is the 

angular frequencies of the ultrasonic wave. Arrows in 
Fig. 3 indicate the temperature being the resonant con- 
dition of WT = 1. The anomaly of (Cn — Ci2)/2 around 
30 K is also well described by the Debye dispersion of the 
solid and dashed lines in Fig. 1. 

In inset of Fig. 3 the temperature dependence of 
the relaxation time r obtained by Cl is presented to- 
gether with the results of (Cn — Ci2)/2 of Fig. 1 
and Cii, that is not presented here. The relaxation 
time due to the rattling motion obeys the tempera- 
ture dependence of r = TQexp{E/k^T) with an at- 
tempt time To ~ 8.8 x 10^^^ sec and an activation en- 
ergy E = 168 K. Utilizing a harmonic oscillation of 
C,{z) = (I/ttzo)^^^ exp(— 2^/2^0), we estimated a mean 
square displacement zq — (l/27r)(/iro/M)^/^ — 0.079 nm 
for Pr-ion in the present potential of the cage^\ This 
result is comparable to zq — 0.048 nm of Ce3Pd2oGegii^ 
and zo = 0.012 nm of La3Pd2oGeeii. 

The twelve Sb-atoms have a distance 0.3542 nm from 
the center of the cage. Because the Sb-atom is absent 
along the [100] axis, the Pr-ion may favor an off-center 
motion over six minimum points of potential at ri = 
(a, 0,0), r2 = (-a, 0,0), = (0,a,0), n = (0,-a,0), 
rs = (0,0, a), rg — (0,0, —a). Here, the mean square 
displacement of Pr-atom extends over the potential min- 



imums as a ~ zo/2 = 0.04 nm. On the other hand, the 
Os-atom locating at 0.4028 nm from the center of the 
cage along the three- fold [111] axis prevents the off-center 
motion along the [111] axis. The Sb-atoms closely locat- 
ing to the two-fold [110] axis may hinder the off-center 
motion along the [110] axis. 

The ultrasonic dispersion has been observed in the 
transverse wave of (Cn — Ci2)/2 of the strain Ev with 
r23 symmetry. The longitudinal modes of Cn and Cl 
consisting of the strain — {2ezz — ^xx — £yy)l^f^ with 
r23 symmetry in part also show ultrasonic dispersion. On 

the other hand, the C44 mode responsible for the elastic 

('2') 

strain with symmetry does not show the dispersion 
effect. These results indicate that the thermally activated 
rattling motion being coupled to the elastic strains Su and 
£ru has the F23 symmetry. It is of particular importance 
to project out the off-center mode for the irreducible rep- 
resentations at a center of the cage with the point group 
symmetry TifS^iS^iS^. Operating a symmetry element R 
of Th on the atomic density pi = p{pi) at the minimum 
point pi{i = 1, 2, • • •, 6), one obtains six dimensional rep- 
resentation matrices Dij{R). The character x{R) being a 
trace of the representation matrix reduces to direct sum 

(2) 

Fi ® F23 © F4 • The projection operator is used to pick 

(2) 

up the Fi, F23 and F4 representations consisting of the 
fractional atomic density of Pr-ion over the six minimum 
points. In the present case of PrOs4Sbi2, the F23 off- 
center mode of ppssu — 2/^5 + 2/96 — Pi — P2 — P3 — P4 
and pr2av = Pi + P2 — P3 ~ P4 with the fractional atomic 
distribution in Fig. 4 is the ground state of the sys- 
tem. pr23u means the distribution of fraction 1/2 at ps 
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FIG. 4: The rattling mode of praa,, = 2p5 + 2p6 — pi — 
p2 — Pi — pi and pr23„ = pi + P2 ~ ps — p4 being responsible 
for the ultrasonic dispersion in PrOs4Sbi2. 

and pe sites and null at pi,p2,P3,P4- And pr23„ is re- 
sponsible for fraction 1/3 at pi and p2, 1/6 a-t Ps and 
Pe, and null at pa and p4. The total symmetric mode 
Pri = Pi + P2 + P3 + P4 + P5 + P6 and the polar mode 
PTi^ = Pi - P2, Pr4„ Pa - P4, Pr4. = P5 - Pe may cor- 
respond to the excited states, ppi represents the mean 
fraction 1/6 over the six sites, pr^^, for instance, has 
fraction 1/3 at pi and null at p2, and fraction 1/6 at 

P3,P4,P5,P6- 

Recent ultrasonic measurements on the clathrate com- 
pound La3Pd2oGe6 by our group has successfully showed 
the rattling and tunneling motions of an off-center La 
ion in cage^''. In the present clathrate compound of 
PrOs4Sbi2 with a cage of Sb-icosahedron, the thermally 
activated rattling motion over the potential hill brings 
about the ultrasonic dispersion with the relaxation time 
T in inset of Fig. 3. With lowering temperature, the 
thermally activated rattling dies out completely with- 
out showing the structural transition. Consequently, the 
tunneling motion of Pr-ion through the hill in keeping 



the site symmetry to be Th is relevant at low temper- 
atures. The tunneling motion being accompanied by 
charge fluctuation interacts with conduction electrons 
through the channels of Pr-Sb bonding. It is remarkable 
that the off-center tunneling motion with r23 symmetry 
may bring about the enhancement of the elastic softening 
in (Cii — Ci2)/2 just above Tc- The theoretical work for 
the interplay of the tunneling motion to the supercon- 
ductivity by Cox and Zawadowski^^ may be relevant for 
the present PrOs4Sbi2. 

In conclusion we have successfully observed the elas- 
tic softening in (Cn — Ci2)/2 and C44 above Tc- It is, 
however, still difficult to determine a CEF state in the 
alternative model of Fzs-rf ^ or Fi-Pf ^ . The field de- 
pendence of the elastic constant in particular is neces- 
sary to settle the CEF state of the system. Nevertheless, 
it is worthwhile to emphasize the fact that the soften- 
ing of (Cn — Ci2)/2 and C44 indicate a crucial role of 
the quadrupolar fluctuation to the heavy Fermion super- 
conductivity in PrOs4Sbi2. Furthermore, the ultrasonic 
dispersion due to the r23 rattling motion with activation 
energy E = 168 K has been foimd. The r23-type charge 
fluctuation associated with the off-center tunneling mo- 
tion in particular may enhance the elastic softening of 
(Cn — Ci2)/2 just above Tc. The more accurate investi- 
gation is necessary to clarify the CEF state in PrOs4Sbi2 
and the interplay of the quadrupole fluctuation and the 
off-center motion of Pr-ion to the unconventional super- 
conductivity. 
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